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A STUDY OF BLAST LO.\DIKG 



TRANSMITTED TO BUILDING FRAMES 

I. INTRODUCTION 
A. Intrcxiuctory Statement 

The loading which a structure experiences due to the pressure waves accom- 
panying blasts frcm high explosive and atomic bcmbs has been the subject of 
much study and research since the use of such bombs was begun In the last war. 
The advent of more powerful atomic bombs and of the "super" bombs has Increased 
even more the urgency of discovering with greater exactness the nature of blast 
loads on structures. 

A recent study In the field of blast loadlng^^^ has concluded that a gener- 
alized triangular pulse Is a good approximation to the actual transient loads 
which the front wall of a structure experiences. However depending on the 
action of the wall covering, the shape of the load- time curves that reach the 
building frame Itself may be quite different frcm the pulse which Impinges on 
the w£lL 1 panels. It Is the purpose of this thesis to study the load-time rela- 
tionship for the load transmitted by the wall covering to the supporting frame. 

The problem is necessarily complex and Is approached by making use of cer- 
tain simplifying assmptions. The most important are: 

(l) Given a typical bay of a one- story frame- type building, consider that 
the wall covering of the bay consists of a "nest" of closely- spaced 
beams of unit width which are simply supported and act Independently 
of one another. (Figure l) These beams are used in this study to 
approximate the action of a slab covering. 

^^^N. M. Newmark, "An Engineering Approach to Blast Resistant Design", 
Proceedings ASCE, Vol. 79 , Sepprate No. 306, October 1953, 
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(2) Assume that the pressure loading which impinges on the wall panels 
can be represented by a triangular pulse function. (Figure 2) 

(3) Assume ideal elas to -plastic action of the beams. (Figure 3) 

(1*) Assxane that the frame members supporting the wall covering remain 
fixed in space^ at least d\iring the interval for which this study 
is made. 

It is doubtful that movement of the frame would have much effect on the 
loads reaching the frame members. A deflection of the frame away from the 
blast source would result in a slight decrease in the beam end reactions; 
but, since the period of vibration of the supports is much larger than that 
of the wall iwinels, the effect shovild be negligible except possibly for load 
pxilses of very long dvtration. 

If the beams which approximate the wall covering were perfectly rigid and 
massless, the shai)e of the load transmitted to the frame would be identical 
to that of the blast pressure pulse (Figure 2). However in our case, the beams 
have mass and also have elas to -plastic characteristics; therefore, the load 

which reaches the supjxjrting members is much altered. 

/ 

In order to determine the frame loeidlng for the given assumptions, an 
elastic analysis and a simple plastic analysis are used. The elastic analy- 
sis is derived from the familiar dynamic beam equation (l) and the plastic 
one comes from a simplified consideration of beam action when a plastic hinge 
forms at the center. From these two analyses expressions for the dynamic end 
reaction of the beams on the supports are derived. 

Frame load charts giving the variation of dynamic end reaction (frame loaui 
per unit height) with respect to time are drawn for different beam resistance 

M. Newmark, "An Engineering Approach to Blast Resistant Design", 
Proceedings ASCE, Vol. 79> Separate No. 306, October 1953 • 
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characteristics euid different load pulse durations. The frame loads sure com- 
pared graphically with the load which the beam receives as a result of the 
blast pressure wave. Load -time curves for the supporting frame are shown 
only for the dvtration of the loading pulse; however, a method is indicated 
hy which the relation may be found for times beyond the end of the pulse. 

As brought out in a subsequent section of this study, the elastic reaction 
expression involves the sum of a trigonometric series. For certain durations 
of loading pulses, the curves for one term, three terms emd six terms are 
drawn for comparison. By means of this graphical comparison of a few repre- 
sentative pulse durations, it is demonstrated that the use of only the first 
term gives resiilts which are fairly accurate for long pulses but not very ac- 
cvirate for short pulses. For all the remaining pulse durations, only the 
first mode (one term) is used although this is only an approximation to the 
exact relation which theoretically would require an infinite ntimber of terms. 
No attempt is made by analytical means to evaluate the accviracy obtainable for 
various pulse dxjratlons by using only the first mode. Such an evaluation 
could be made the subject of fxirther study. The eidvantages of using the first 
mode approximation in the construction of the frame load charts is apparent, 
however, emd such use is intended. 

B. Notation 

Terms used in the following sections are defined as they appear. They 
are assembled here for convenience in reference. Dimensions of the quantities 
are indicated by the symbols F, M, L, T in parentheses meaning respectively 
force, mass, length and time, 
a^ = El/m 

9 = angular displacement 

9 = angular velocity (1/T) 
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T 

n 

V 



angular acceleration (i/T ) 



=r 



pdt (FT/l) 



Impxilse of loeid. 
time distribution factor 

l/(ti/Ti) 

length of beam (wall panel strip) (l) 

mass per unit length of beam (m) 

dynamic bending moment at center of beam (FL) 

fully plastic moment of beam (FL) 

maxiimmi static center moment = /8 (FL) 

ultimate moment - carrying capacity of beam (FL) 

instantaneous value of applied pressui’e load (f/l) 

m 

= space distribution factor 
maxitman value of applied pressure load (F/l) 
dynamic end reaction (f) 
maximum static end reaction = P X/2 (F) 

time (t) 

duration of load p\ilse (T) 

time of beginning of second elastic phase (T) 
time of beginning of plastic phase (t) 



t-tp (T) 

(T) 

natureil periods of vibration of beam (T) 
fundamental natural period of vibration of beam (T) 
a function of time alone 
shear (F) 

deflection of beam (l) 



w 
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w a linear velocity of beam (L/T) 

V = distance from center of mass to axis of rotation (l) 

V r= acceleration of mass center (l/T ) 

Wg = effective center deflection of beam at onset of second 
e 

elastic phase (l) 

w * maximum center deflection of beam which occurs at t = t (L) 
c e ' 

m 

Wy a center deflection of beam at t = t^ (l) 

X a characteristic function for beam 
n 

= natural (circular) frequencies of beam « 

= fundamental natural frequency of beam (i/T) 




II. INITIAL ELASTIC BEHAVIOR OF WALL PANEIS 



A. Statement of the Problem 




When a beam resists moving loads ; the governing differential equation of 
motion for the beam is, according to the ordlnsiry Bernoxilli -Euler theory of 
flexure, 

' mm 
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a « p(^) 
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y x„ 0(t) 
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m 



(1) 



where a a El/m; w a deflection measured normal to the beam axis; m = mass per 
unit length; q » external applied unit load; x » distance measured from origin 

at left end of beam; and t = time. The external applied load is defined by a 

m 

space distribution factor, p(x) * z ^n ^n ^ time distribution factor, 0(t). 

n 

^^^See for example, S. Timoshenko, "Vibration Problems in Engineering", D. Van 
Nostrand Co*, Inc«, New York, N. Y., Second Edition, July 1937; P* 332* 



Figure 4 is a simplified sketch of the loading on a typical beam. For the 
given triangular loading pulse, p(x) = 1 since the beam is assumed to be uni- 
formly loaded, and )Z5(t) = P^(l - t/t^). The characteristic function, X^, for 
a simply supported beam is sin (nnx/^). The constants, P^, are defined as: 



P 

n 




dx 



( 2 ) 



Substituting in (2) and performing the Integrations leads to: 

^ 3' 5, 7 

P = 0 for n = 0, 2, 4, 6 

The solution to (l) is taken in the form 




n 



(3) 



where X = sin (nnx/£ ), a f vine t ion of x alone, and where x is a function of 
time alone. The substitution of (3) in (l) with the subsequent separation of 
vsiriables leads to the following differential equations : 



1 t I t 



X 
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(^) 



• • 2 

T + U) X 

n n n 



Pq 0(t) 

m 



(5) 



The solutions of (4) and (5) are respectively: 



7 

= sin (always true for simply- (6) 

supported beam ) 



T *5 A cos u t + B sin u t 
n n n n n 



n 

mu) 



sin 



n 



(jj t 

n 



t 



0(t) 



COS 



tj tdt 
n 



cos 



“n‘/ 



0(t) sin 



u tdt 
n 



(7) 



Equation (7) is the solution to the single-degree of freedom Equation 
The first two terms in (7 ) axe the free vibration terms and the last two the 
forced vibration terms. For initial conditions of zero displacement and zero 
velocity, the free vibration terms disappear since 
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dx 



B 
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I “ o <*) =<0 

2 

I 

»o 



dx 



where v^(x) and w^(x) are the displacement and velocity respectively at zero 
time. 
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S. Timoshenko, 
New York, N. Y 



"Vibration Problems in Engineering", 
,, Second Edition, July 1937, P« 3^* 



D. Van Nostrand Co., Inc., 
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Substituting Equations (6) and (7) into the expression for deflection (3), 
yields the following: 



. n . niix 

“ ‘ ^ T 



-t 

sin / 0(t) cos u^tdt - cos T 0(t) sin w tdt 

Jo Jo 



(8) 



Substituting the values of and 0 (t) in (8) eind integrating, results in the 
following expression for w: 



w = 



4 1 m . nnx 



nn mu u 
n n 



sxn-^ 



(1 - cos u t) - r— ( 1 

“hi 



n 



sin u) t 



u t 
n 



2 2 2 ■ 
Noting that u^ for a simply- supported beam = (n n/X ) ^, vhere a =^El/m: 



w = 



m 

tT^EI 



1 . nTTx 

,5 X 



n=l,3,5--- 



(1 - cos u t) - ~ I : 

° h \ 



sin u t 
n_ 

u t 
n 



(9) 



B. Reaction eind Moment Expressions 
1. I)ynemlc End Reaction 

The shear at any point x along the beam is given by: 



V = - El ^ 



Performing the differentiation on (9) and evaluating the expression at the 
end of the beam (x = O), the dynamic end reaction is: 



9 







Sin (J t 
n 



(1 - cos u^t) - ^ I 1 - 

^ \ tJ X 

1 \ n 



If the maximum value of the external load, were applied as a static load, 
the maximum static end reaction woxild he: 



P 1 

R - P 
■ 2 



For ease in plotting, the reaction expression is put Into dimensionless form 

2 2 
by noting that = 2nn /T^ since = T^^/n . 



00 




n«=l,3>5“ 



- l/(t^/T^) • t/T 




sin 2"n^ t/T^ 
2iin^ t/T^ ^ 



( 10 ) 



2. Magnitude of Reaction for Long Pulses 

As the dxiration of the load pulse goes to infinity, the coefficient of the 
second term in (10 ) approaches zero and the reaction approaches: 



00 




n.1,3,5 — 



VTi 



The maximum reaction then approaches: 



10 




The limit of the summation exists as shown. It is given by the Riemann Zeta 
Function, the values of which are tabulated. ^ 

For X o 2, ^ (2) = n^/6 = 1.645; or: 



t,/Ti 




16 

2 



n=odd 



(1 - lA) T 



2.00 



Then for pulses of duration approaching infinity, the dynamic reaction ap- 
proaches a value twice the maxi mum static reaction. This resvLLt is to be 
expected from energy considerations. 

3. Magnitude of Reaction for Short Pulses and Constant Imptilse 

For short pxilses of duration approaching zero, the reaction should also 
approach a limit. This limit can be fouind from a consideration of the reac- 
tion expression which holds for time beyond the end of the pulse as the pulse 
dixration becomes very small and total Impulse remains constant. Figure 5 
shows the new time origin at t s= t^. 

As before, the deflection, w = ^ X t: . However for t > t, , the forced 
^ ^ n n 1 

n 

vibration terms in Tq drop out, leaving: 

T = A cos (j t + B sin u t 
n n 1 n n n 

Jahnke and F. Elude, "Tables of Functions with Formulae and Cxurves", 
Dover Publications, New York, N. Y., Fourth Edition, 1943* 
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The initial conditions of displacanent and velocity for t = 0 can be fovmd 
from Equation (9)- Setting t = t^ in (9): 



kP I 
m 



/ \ m \ 1 . nnx 

n=odd 



sin u t, 
n 1 

u t, 
n 1 



cos u t, 
n 1 



Differentiating (9) with respect to time and setting t = t^: 
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4p i' 
m 



IT^EI / n 

n=odd 



i, slniHi 
5 I 



COS u t, , 

. . nil 

u sin u t, + - — 

n n 1 t^ 



Values of A and B are then: 
n n 
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Jo "o**' \ 



dx 
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Jo n 



, ^ . nnx 

where X = sin —r— 



dx 



A = 




sin (J t, - u t^ cos (J t, 
n 1 n 1 n 1 

w t^ 
n 1 



1 Jo 
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(x) X dx 
o' n 
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n 1 



cos 



u t- 
n 1 



u t, - 
n 1 



Then for t s t - t 
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sin 
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u t 
n 



u t, - 
n 1 



ui 
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n 1 
n^l 



cos 




+ sin t 
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/ u) t, sin u) t, + cos (j t, 
[ n 1 n 1 n 1 

\ n 1 




As before, find the shear and evalxmte at x *> 0 in order to get the dynamic end 
reaction; 



CO 




For any pulse duration, the impulse of the load is: 
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■^1 

pdt . 1/2 
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Now let t, — >■ 0 keeping the total impulse constant. P then increases 

1 ID 

rapidly, and in terms of I is given by: 



P 



m 



21 



t 



1 



Replacing P^ in (ll) by its value in terms of I and taking t^^ Inside the sum 



mat ion: 




n=odd 



* 

cos u t 
n 



sin u t^ - u t, cos u t, 
n 1 n 1 n 1 



u t, 
n 1 



+ sin tJ t 
n 



„ / t, sin u t- + cos w t- - 1 

* / n 1 11 1 n 1 

u t ^ 

n 1 



Now let tj^ approach zero: 



lim 

0 



R _ ^ 



1 “n • klJL 

=-; 2 - 

n Ti 



n * 

^ sin u t 
2 n 

n 



n=odd 



n=odd 



Noting again that = 2n/T^ and replacing I by (l/2) P^ t^^; 



Bg- r 

0 



sin u t 
n 



n^ 



n=odd 



Take the first mode by letting n = 1: 



14 



t, — >0 



^ “ I ^ 2" I;; 



( 12 ) 



The foregoing derivation has shown that for the conditions of very short 
pulse duration (say when ^ l/ 3 ) smd constemt impulse, the TnavlTmim elastic 

reaction approaches a value equal to (8/ it) A plot of reaction for 

t > tj^ in which the load pulse satisfies the two conditions shovild verify ( 12 ). 
4. Dynamic Center Moment and Beam Resistance Parameter 
The mcsnent at any point x along the beam is given by; 

^2 
o w 

M » - El 



Performing the differentiation on (9)» evaluating at x = ^/2 and noting that 
2 

P / /8 is the maximum static center mcment if P were applied as a static 

m 

load: 

y -— ■ n-1 

Lii " 

^ TT^ j n^ 

n=odd 

- • t/T^ 

A plot of Mjj/Mg versus the time parameter t/T^ furnishes a convenient means 
of defining the limits of elastic and plastic beam action for different pulse 
durations. As in the reaction plots, the first mode (n = l) in ( 13 ) is used 
as a good approximation. Given a cujrve of M^/Mg versus t/T^, the fully plastic 



1 - cos 2 nn^ ~ 



1 - 



sin 2 TTn^ t/T^ 



2TTn t/Tj^ 



( 13 ) 
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mcanent of the beam cross-section may be considered to have been reached for 
different values of M^. When a given Is said to equal Mp/Mg, the time 

corresponding to the particular value of is the time of hinge formation, 

tp. The beam resistance parameter, Mp/Wg “a-y thotjght of as the ratio of 
the yield value of the beam's moment reslstaince to the majclmum mcment caused 
by the static application of P^. 



Ill PLASTIC BEHAVIOR OF WALL PANELS 
A. Statement of the Problem 
1. Assumptions in the Simplified Plastic Analysis 

In order to derive a fairly simple expression for beam reaction when the 
beam ceases to behave elastically, the following assumptions are made: 

a. A plastic hinge forms at the center of the beam when the dynamic center 
moment, reaches a value eqxial to the fully plastic moment of the sec- 
tion, Mp. This time is designated t^. 

b. Neglect the elastic bending of the beam. 

c. Assiune that deflections are small. 

d. Assione that remains equal to Mp until the beam velocity becomes 
equal to zero. At this point, becomes less than Mp, the center 
hinge disappears, and the beam behaves elastically again. 

FMgure 6 represents the idesilized plastic behavior of the beam conforming 
to the above assumptions. 

Figure 7 is a free body diagram of half of the beam from which equations 
of motion may be written as follows: 




where w 
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(14) 
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where 



1 m/ fiV m£^ 
3 \2 J “ 




^ M_ 

TT ■ “p 



”2r 



(15) 



B. Reaction Expression 

Frca Eqviatlon (15) solve for the angular acceleration: 




( 16 ) 



Substitute (l6) Into Equation (l4) and solve for R: 



R 






where 



p » ftuactlon of time. 



(17) 



In Equation (17)> P is the external load acting on the beam at time t' measured 
after the formation of the center hinge. Figure 8 represents the time origins 
for the different phases of beam action which occur during the lo£id pulse. 

With 



P = p(t') = P, 



m 




and t’ ® t - t , 
P 



and noting again that Pjjj>^/2 is the maximum static reaction: 



17 




( 18 ) 



Differentiation of (l8) with respect to the time parameter, t ' gives the 
result that the slope, S, of the plastic reaction curve depends only on the 
pulse duration, or: 

S = - J f if K = l/(t^T^) (19) 



C. Duration of Plastic Behavior 

1. Expression for t^ 

As stated in the assuBiptions for the plastic analysis, the center hinge 
disappears when the beam velocity becomes equal to zero. The time, t^, at 
which the beam stops moving down is therefore the end of the plastic phase 
eind the beginning of the second elastic phase. 

In order to find t , an expression for the angxilar velocity of the beam 
can be derived which holds during the plastic phase (t^ Equating 

this angular velocity to zero and solving for the time will yield the equa- 
tion for t . 

e 

Differentiating Equation (9) with respect to time and evaluating at 

gives the linear velocity of the beam center, w^ during the initial 
elastic phase. For small deflections, the average angular velocity of half 
of the beam wovild be = ZjJi w^. Performing these operations and setting 
t = t^ will give the equation for the initial angular velocity at t' = 0 
(t = tp). (Figure 8) Thus; 
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n-1 r- 



u sin u) t 
n n p 



^ (l - cos u t ) 



n p‘ 



( 20 ) 



The required angular velocity may now be fovind by Integrating (l6) which 
was derived from the equations of motion for the plastic phase. 

Noting that p in (l6) = 



P(f) » P, 



f + t 
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1 - 



O' = 



3P 
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m i 



t' + t 
1 - -T H 






Integrating and evalviating the constamt of integration by using (20): 
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u) sin <J t - ^ 


1 1 - cos u t ) 




n n p t^ 
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Simplify by letting n = 1 (first mode) and note that = ir^El/mX^: 
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3nu^ 



sin 




(1 



cos u), t ) 
1 P 



( 21 ) 



When (21 ) is equated to zero, the resulting equation is a quadratic in t' 



which con be solved by the quadratic 
put in dimensionless form: 




formula giving the following resvilt when 





IV. SECOND PHASE OF ELASTIC WALL PANEL BEHAVIOR 
A. Statement of the Problem 
1. Application of Dynamic Beam Equation 

Equation (22) defines the time at which the beam again behaves elastically. 
The derivation of this second elastic reaction expression follows the same 
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20 



method as used on the Initial one; however, one free vibration term is added. 
When t" = 0 (t = tg), the initial velocity, v^" of the beam is zero, but the 
Initial displacement, w^" has a finite value. An examination of the formulae 
for and (Section II ) shows that the first free vibration term in Equation 
( 7 ) has a value other than zero. 

For time measured from t (t" » t - t ), the elastic deflection is still 

e e 

of the form: 

«• 

W"(x)= 

n 



where = sin niix/Z and 
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A cos 
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U) t” 
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mu 
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cos w t’*dt** 
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- cos 






(23) 



where = U/nn when n = odd. (See Equation (2).) 

The constant A is defined as 
n 
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Jo 



(x) X dx 
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x/dx 
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f w (x) X 
/ o ^ ^ n 

wo 



dx 



(2U) 
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Now let w "(x) = w sin -7- where w is the effective center deflection at 
o c / c 

e e 

the beginning of the second elastic phase. At first glance it appeetrs that a 
good approximation of w^"(x) should be given by 



w ”(x) c w sin ~ 
o c p 

m ^ 



where w is the center deflection of the beam in the plastic phase at the 
instant the beam reaches its maximum deflection in the positive direction (in 
direction of applied load) and begins to reboxmd elastically. However^ it is 
apparent that the maximum amplitude of the elastic center deflection is much 
less than the corresponding maximum amplitude of the plastic center deflection. 
Figure 9 illustrates roughly this difference. The solid curve gives an idea 
of the actual variation of reaction with center deflection throughout the two 
elastic phases and the single plastic phase. The dotted curve is an approxima- 
tion to the act\mil behavior which suggests a simple means of obtaining a close 

estimate of w .As shown in the figure, the maximum plastic center deflec- 
'^e _ 

tion, w is much greater than w ; and v , the center deflection at the time 
c c y 

m ^ _ 

of hinge formation (t ) is approximately equal to w 

^ e 

The final expression and best approximation to be used for w^"(x) is thus: 



nx 



w "(x) = w sin ^ 
o'-' y / 



The value of w may be fotind from Equation (9) by setting x = //2 and 

y 

t s= t . The result is: 
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(1 - cos u^tp) 




(25) 



n=odd 
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Nov in X 

n 



sin — (24), let n = 1. 



A then reduces to: 
n 



A 



n 




. 2 TTX , 

V sin -j- dx = V 

y X y 



The final expression for A^ is thus identical to and is given by (25)- 

Go back to the expression for Equation (23), and substitute the value 
of 0(t") 



0(t") = P 



t" + t 



m 



After carrying out the integrations in (23), combining it vith and simpli- 
fying, the expression for elastic deflection for < t < t^^, is: 



v"(x) = 




nTTX 

T 



cos (j t" 
n 




cos u t ) 
n p" 




sin (J t 
P P 

u t 
n p 





B. Reaction Expression 

Having derived Equation (26) for deflection, the reaction expression for 
the second elastic phase may be found as before by evaluating the shear at 
X = 0. In dimensionless form for plotting, and using the first mode only, the 
reaction equation for t < t < t , is: 
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sin 2n t"/T^ 




(27) 



V. RESULTS AND DISCUSSION 
A. Frame Load Charts 

1. Construction of Charts 

In the previous sections of this paper, expressions for the beam end re- 
actions have been derived for the elastic and plastic phases of beam behavior. 
The reaction is expressed as a ratio of the dynamic end reaction to the maxl- 
m\jm static end reaction, (Rj^/Rg); the time parameter is written as the ratio 
of the elapsed time to the fundamental natural period of the beam (t/T^); and 
the beam resistance parameter is expressed as the ratio of the fully plastic 
moment of the section to the maximxmi static center moment (Mp/Mg). For con- 
venience, the following table lists the various phases of beam action, the 
time range of each phase (Figure 8), and the governing reaction equation: 



Phase 


Time Range 


Equation 


Initial Elastic 


o < t < t 

P 


(10) 


Plastic 


t < t < t 


( 18 ) 




P e 


Second Elastic 


t < t < t, 
e 1 


(27) 


Figures 23, 27, and 


29 contain plots of all three phases. 


The remaining 



reaction curves include only the initial elastic and the plastic phases. On 
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the same graph with each reaction plot is shown the external load which im- 
pinges on half the beam. Thus at t = 0, the load on half the beam is P^/^/2; 

the mAyinnnn static end reaction is slIso P z^/2; therefore the initial ratio 

ni 

of the half -beam load to the maximum static reaction is equal to one. The 
external load line thus starts at unity at t = 0, and drops in a straight 
line to zero at t = t^^. This is, then, a graphical representation of the load 
which arrives on half of the beam, and what this half beam transmits on to 
the frame as end reaction. 

The moment cvirves are vised to find the vaJue of t^ for a given pulse 
dvuration t^^, and resistance parameter (Mp/Mg) . (Figure l6) The duration 
of the initial elastic phase is thus determined graphically from the plot of 
Mp/Mg versus t/T^^ by allowing Mp to equal Mp as Mp changes over a range of 
values. This brings in the shape and physical characteristics of the beam 
since the fully plastic moment varies with the shape and material. 

The duration of the plastic phase is fixed by Equation (22). The loci of 
these times as found by solving (22) form a smooth curve and are plotted on the 
reaction graphs. It will be noted that for a given pulse duration, the loci of 
Equation (22) indicate the expected results; that is, that beams of low resist- 
ance (smaUu Mp/Mj,) have long plastic phases and vice versa. 



In the construction of the reaction curves, it was found that a discon- 
tinuity always existed between the elastic eind plastic phases. This "Jump" 
in reeu^tion value is due to the relatively rough assumptions of the plastic 
analysis used, the main error being the neglect of elastic bending of the beam. 
Since the source of the discontinuity was known, the reaction values were cor- 
rected by the amount of the initial reaction error. 

Plots of the second elastic phase were made only for a few representative 
find 



Mp/Mg as follows: 



values of 
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t/Ti 


V«3 


Fig. 


a. As an example of a short pulse 


0.8 


1.0 


23 


b. As examples of long pulses 


5.0 


1.4 


27 




5.0 


1.0 


27 




10.0 


1.4 


29 




10.0 


1.0 


29 



2, Coinpaxlson of First Mode and More Exact Reaction Curves 

As stated in the introductory statement, no attempt is made in this study 
to find out by analytical procedures the accuracy with which the first mode 
approximation represents the exact elastic reaction curve. However, for the 
initial elastic phase. Equation (lO) was expanded and data was collected to 
plot curves for one term (n «= l), three terms (n = 1, 3 and 5)» and 6 terms 
(n = 1, 3# 1 } 9 and 11) in order that a graphical cc*ai>arison could be made. 

Figure 11 shows the results of the plot for a very short pxalse. It is 
obvious that the first mode approximation is not good for load pvilses of such 
short duration. Figure 35 illustrates the very great Improvement in the ac- 
curacy of the first mode approximation for pulses of very long duration. It 
would seem reasonable to expect then, that as the load pulse duration increases 
from a small to a large value, the accuracy of the first mode shotild show a 
gradual improvement. Figure 19 furnishes evidence of this trend since it shows 
the first mode curve to be a better approximation than it is for = 0.1 

but not as good as for = 100.0. 

^^^Thls subsection refers to the unconqpensated first mode reaction curves as 
derived from Equation (10 ) by setting n = 1. These curves are indicated in 
Figures 11, 19 and 35* As brought out in the following subsection, all other 
reaction plots are constinicted using the first mode multiplied by a compen- 
sating factor. 
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One aitxiation in which the first mode curve gives a false impression is 
encountered at very small values to t/T^ and Mp/Wg. An example is Figure 11. 
Here the reaction in the plastic phase for Mp/Mg = .05 appears to become neg- 
ative before the end of the pulse. It is apparent, however, that negative 
values woiild not be indicated if the plastic phase plot were started from the 
more exact initial elastic curve (for 3 or 6 terms). An examination of pulses 
of longer duration will indicate this same ambiguity at very small values of 
t/T^ and Mp/Mg. 

3. Method Used to Moke First Mode Approximation More Conservative 

Figures 11, 19 , and 35 show that the first mode reactions plotted from 
Equation (lO) generally furnish values lower than those indicated by the more 
exact curves. In order to retain the advantages of the first mode approxima- 
tion but at the same time to obtain more conservative end reaction magnitudes, 

the reaction expressions for the two elastic phases (Equations (lO) and (27) 

2 

with n s: l) were multiplied by a factor of n /8. This step also increased the 
plastic reaction values by the same factor since the plastic phase merges with 
the two elastic phases. The reasoning behind the use of this factor follows. 

The use of the first mode curves for elastic reaction neglects the contri- 
bution of the higher modes in (lO) and (2?) to the end reaction. Svich a sim- 
plifying approximation may be considered as being equivalent to using an in- 
stantaneovis beam loading represented by the first term in the Fotirier Series 
expansion of a uniform load. This first term in the Fourier expansion is merely 
a simple sine curve distribution of load. Therefore, the actual pressure load 
which at any instant is uniformly distributed, is replaced by a sine curve load- 
ing the total magnitude of which is very close to that of the actual load. 

Figure 36 illustrates the shapes of the actual and approximately -equlveulent 



loadings 
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The ratio of the end reaction cadculated from the correct loading to that 
derived from the equivalent sine curve loading is found as follovs: (Figure 36) 



"(1) = ^ 



"(2) P*'"' slo ^ dx = 



then 

"(2) ■ S' 

It seems reasonable then, that reaction values more nearly eqvtal to those 

resxilting from a uniform load will result if Equations (lO) and (27) (setting 

2 

n = 1) are multiplied by the factor n /8. Such a procedure was followed in 
this paper. 

Figures 11, 19 and. 35 contain plots of the initial elastic reaction expres- 
sion wherein a comparison is made of the curves plotted from Equation (lO) for 
the following conditions: 

(1) n = 1 

2 

(2) n = 1 using multiplication factor of n /8 

(3) three terms (n = 1, 3, 5) 

(4) six terms (n = 1, 3, 5, 7, 9 , H) 

This graphical comparison shows that the use of the first mode multiplied by 
2 

the factor 11 /8 gives results which are more conservative approximations to the 
more exact curves. 

Since vadues of t are determined graphiceilly from the M^y^Mg versus t/T^ 
plot, another source of error in the first mode reaction curves is the use of 
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only the first term (n » l) in Equation (13 )• However, the elastic center mo- 
ment varies as ~ for n equal to odd integers; therefore, the approximation for 
n^ 

center moment is less in error than the corresponding approximation for reaction 
since elastic reaction varies as ~ . A comparison of the center moments due to 

Tl 

the two loadings in Figure 36 gives the following result: 



M 



( 1 ) 



. 0.125 p(t)/^ 



then 



M 



( 2 ) 



= (.500 - .182) = 0.129 p(t)/‘ 

TT 



M^2) -129 ■ 1.032 



.970 



By the foregoing reasoning, more accxirate values of t /T^ could be gained 
by m\Jltiplying Equation (13) by 0.970. This factor is so close to unity, how- 
ever, that little change in the values of t would result. In this study. 
Equation (l3) with n = 1 was used as derived without applying the correction 
factor. 

4. Discxission of Different Phases of Wall Panel Action 

In the introduction to this paper, the statement was made that the charac- 
teristics of the wall panels would cause the shape of the load transmitted to 
the frame to be much different from that which impinges on the panels themselves. 
A study of the Frame Load Charts (end reaction curves) bears out this statement. 
The wall panels are acted upon by an initieil-peali, triangular loading pulse 
but send to the frame by means of the end reactions, a loading which has been 
radically altered in shape. 
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The unit load which reaches the supporting members rises gradually from 
zero along the initial elastic reaction curve. Depending on the duration of 
the blast pulse and the beam resistance parameter, the beam (wall panel) be- 
comes plastic at a certain time emd the reaction then decreases linearly with 
time at a rate depending only on the pxilse duration. Then if the panel has 
sufficient elastic reslsteince, the action eigain becomes elastic and remains so 
throughout the remainder of the load pulse. 

For very short pulse durations, the plastic phase continues to the end of 
the pulse and no second elastic phase occurs. Figures 11, 13, and 15 illustrate 
this behavior. As the pulse durations increase, the second elastic phase be- 
gins to appear starting with panels of high elastic resistance (large Mp/Mg). 
Figures 2? and 29 show how (for long pulse durations) the two elastic reaction 
curves oscillate about the Incident load line. It should be noted here that 
damping will tend to decrease the amplitude of the elastic oscillations about 
the load line. If the effects of damping could have been included in the de- 
rivations of the elastic reaction expressions. Figures 27 and 29 would proba- 
bly indicate a gradual decrease with time in the amplitude of the vibrations of 
elastic reaction about the extenml load line. Thus, in the second elastic 
phase for low beam resistances, the elastic reaction would become coincident 
with the load line as t approaches t^^. 

Theoretically, if a wall panel had Infinite elastic resistance, its end re- 
action would follow the initial elastic curve and the frame loads would rise 
to a maximum. Plastic action telces place however, since the panels have 
finite elastic resistance, and prevents the frame loads from building up to 
these peak values. The formation of the center plastic hinge, thus effectively 
cuts off the top of the elastic reaction curve and consequently permits sma ller 



loads to reach the frame. 



I 
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The occurrence of negative values of reaction in the elastic oscillations 
about the load line may be highly significant under certain conditions. If 
the fasteners tised to attach the vail panels to the supporting frame are not 
designed to withstand tension loads, it is very possible that such fasteners 
will fail at the time the wall panel end reaction first becomes negative. 

Such failvure would alter the frame loads in the case of load pulses of suf- 
ficient duration > 0.8 in the frame load charts), limiting the loads 

transmitted by the panels to the supporting members to a single positive pulse. 



B. Tables and Graphs of Results 

Table 1 is presented for convenient reference. For given values of pulse 
duration (t^/T^^), plots are adso made of the beam resistance i>etrameter (Mp/Mg) 
versus both the duration of the initial elastic phase the duration 

of the plastic phase (t^/T^ - t (Figures 37 and 39) Similar graphs are 
drawn with Mp/Wg as the fixed parameter. (Figures 38 and ho) 

Figure 37 shows that the d\iration of the initial elastic phase increases 
almost linearly with increasing values of Mp/Mg for fixed load durations. 

Figure 38 demonstrates that for a given beam resistance, the duration of the 
load pulse has little effect on the length of time which the panel remains 
elastic except in the case of the shorter load pulses. 

As might be exi>ected, there is a very great increase in the dxiratlon of the 
plastic phase (for given pulse duration) with a decrease in wall panel resist- 
ance. (Figure 39) The increase is most pronovinced for long load pxilse dtira- 
tions. 

Figure 40 shows that the effect of load duration is largest for low values 



of x>anel resistance, and that, in general, the plastic phase duration decreases 
as the duration of the load decreases (with given panel resistance). 



Table 
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C. Frame Loads After Termination of Load Pulse 
For time Toeyond the end of the pvilse, the derivation of the elastic reaction 
expression is similar in procedure to that used in Section II in finding the 
limiting reaction for short pulses and constant load impulse. Eqxiation (26) 
wo\ild he used to find the displacement and velocity at t » t^^ in order that 
A_ and B could be evaluated; the expression for w would then be found where 



w a 




a 

n 



and still equals sin vm-xj £, . Reaction could be found as before by evaluating 
the shear at x = 0. 



VI. SUMMARI AND CONCLUSIONS 
A. Applications of Frame Load Charts 

A method has been fo\ind for determining the load -time relation for blast 
load transmitted by the wall panels of a stnicture to the sujjporting members. 

The frame load charts may be used for three possible types of wall panel be- 
havior; a purely elastic action, an elasto -plastic action, and a brittle action. 

For a given pulse duration, and purely elastic behavior, the initial elas- 
tic phase reaction curves apply throughout the duration of the pvilse. This 
assumes, of course, that the wall panel has sufficient elastic resistance (high 
Mp/Mg) to behave in this manner. 

A more usual wall panel behavior is the elasto-plastic type. Here the elas- 
tic resistance is not high enough to confine the action to a purely elastic one. 
Given the pulse dviration and beam resistance parameter, the frame load charts 
used in con^vmction with the dynamic center moment curves indicate the loads 
which are transmitted by the wall panels to the frame. 
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Wall pemelB constrxocted of brittle material will behave elastically xip to 



the ultimate mcment carrying capacity of the panel (m^) amd then fractxire sud- 
denly. The time of fracture could be determined from the dynamic center moment 
curves by setting = My over a range of vedues of My. The frame loewi charts 
would indicate a rise in load along the initial elastic curve xip to the point 
of fracture at which time the reaction would drop to zero. 

In order to xise the frame load charts for all three types of wall panel be- 
havior, the following information must be available: 

a. Wall panel span aiid a curve of side -on overpressxxre versxis time 
from which P^, t^, Rg, and Mgj can be calcxilated. 

b. Dimensions and properties of wall panel in order that T^, Mp 
and My nay be calcxxLated or estimated. 

Entering the charts with and Mp/^^ or My/Mg, the magnitxide and time 

dlstribxrtion of the xiniform frame load can be determined. 



B. CONCLUSIONS 

The frame load charts presented in this stxidy represent an attempt to eval- 
uate the loads received by the supporting members from the building covering 
when the walls are subjected to a blast pressxnre wave. Rough assxm^jtions had 
to be made in some cases in order to obtain fairly simple e:Q)ressions from which 
the load charts coxild be constrxmted. Also, since the elastic reaction expres- 
sions involved the sxm of an infinite series of terms, they could not be used 
in their exact form. Therefore, the first mode approximations to the exact 
elastic reaction expressions were modified by a compensating factor of ti /8 and 
then xised to construct the frame load charts. It is felt that this step, 

M. Newmark, "An Engineering Approach to Blast Resistant Design", 
Proceedings ASCE, Vol. 79> Separate No. 'i06, October 1953* 



eilthovigh not theoreticaJLly exact, provides good engineering approximations to 
the actual "blast loads which a frame would experience. Also, keeping in mind 
the additional condition that uncertainties in loading and in the structural 
properties of the wall panels will always exist, the author feels that the 
frame loads indicated "by the charts are sufficiently accurate to "be used in the 
design of building frames. 

The basic relations for panel end reaction in all phases of elastic and 
plastic action which occur diuring the loeui pulse have "been derived. For the 
values of piilse duration and beam resistance parameters considered, the load 
charts are complete throvigh the initial elastic and plastic phases of wall 
panel action, and representative curves for the second elastic phase have been 
drawn. Thus the main pui^pose of the study -- that of determining the shape of 
the loads transmitted to the frame — has been accomplished. Additional work 
is required to finish the second elastic phase and to extend the study to time 
beyond the end of the pulse using the methods outlined in this paper. 
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Figure 1- lypical Bay with Nest of Beams Approximating Wall 
Covering 




Time 



Figure 2- Triangular Load Pulse 
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Figure 3— IdesLlized Elasto-Plastlc Moment -Curvature Relation 



p(t) = (1 - t/t^) 
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Figure 4 — Typical Beam with Loading 



39 




Figure 5- New Time Origin for t > 




Figure 6 — Plastic Behavior of Beam 



- 




i 




ko 




Figure 7 - Free Body Diagram of Plastic Beam 




Figure 8 - Time Origins for Elastic and Plastic Phases 
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Figure 9 - Graphical Representation of Actvial and Approximate Variation 
of End Reaction with Center Deflection 
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Figure 10 - Ratio of Dynamic Center Moment to Maximum Static Center 
Moment Versxis Time Parameter " 0,1) 



m 

I 






I 



I 

I 

I 

I 

I 

I 




**3 




Figure 11 - Ratio of Dynamic End Reaction to Maximum Static End Reaction Veraus Time Parameter 
(t^T^ - 0.1) 
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Figure 12 - Ratio of Dynamic Center Moment to Maxlmm Static Center 
Moment Versus Time Parameter = 0.2) 
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Figure l4 - Ratio of Dynamic Center Moment to Maximum 
Static Center Moment Versxia Time Parameter 

(t^/T^ - 0.3) 
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Figure 15 - Ratio of Dynamic End Reaction to Maximum Static End Reaction Versus Time Parameter 

(t,/T, =0.3) 
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Figure 16 - Ratio of Dynamic Center Moment to Maximum Static Center 
Moment Versus Time Parameter = 0 . 4 ) 
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Figure 17 - Ratio of dynamic TMd Reaction to Maximum Static End Reaction Versus Time Paxameter (t^/T. 
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Figure l8 - Ratio of Dynamic Center Moment to Maximimi Static Center Moment Versus 
Time Parameter = O.5) 
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Figure 19 - Ratio of Dynamic End Reaction to Maximiza Static End Reaction Versus Time Parameter = 0*5) 
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Figure 20 - Ratio of Dynamic Center Moment to Maximum Static Center 
Moment Versus Time Parameter « 0.6) 
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Figure 21 - Ratio of Dynamic End Reaction to Maximum Static End Reaction Versxis 
Time Parameter ■ 0.6) 




Figure 22 - Ratio of Dynamic Center Moment to Maximum Static Center Moment Versus 
Time Pairameter = 0.8) 
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Figure 23 - Ratio of E^maniic End Reaction to Maximum Static End Reaction Versus Time Parameter (t,/T, * 0«8) 
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Figure 25 - Ratio of Dynaanlc End Reaction to Maximum Static End Reaction Versus Time Parameter ■ l.O) 
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Figure 26 - Ratio of Dynamic Center Moment to Maximum Static Center Moment Versus Time Parameter 
(tj/Ti =5.0) 
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Figure 27 - Ratio of Dynamic End Reaction to Maxtnrum Static End Reaction Versus Time Parameter (t^/T^ «= 5.0) 
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Figure 28 - Ratio of Dynamic Center Moment to Maximum Static Center Moment Versus 
Time Parameter “ 10 *0) 
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Figure 29 - Ratio of Dyiiamic End Reaction to Maximum Static End Reaction Versus Time Parameter * 10,0) 
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Figure 31 - Ratio of Dynamic End Reaction to Maximum Static End Reaction Versus 
Time Parameter (t^/T. = 20.0) 
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Figure 32 - Ratio of Dynamic Center Moment to Maximum Static Center Moment Versus 
Time Parameter = ^O.O) 
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33 ~ Ratio of Dynamic End Reaction to Maximum Static End Reaction Versus 
Time Parameter (t,/T, ■ 50«0) 
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Figure 34 - Ratio of Dynamic Center Moment to Maximum Static Center Moment Versus 
Time Parameter = 100. 0) 
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Figure 35 - Ratio of Umanilc End Reaction to Maximum Static End Reaction Vers\iB Time Parameter (t^/T = 1CX)*0) 
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Figure 36 - Con^axlson of Actual and Equivalent Been Loading 
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Figure 37 - Plot of Duration of Initial Elastic Phase (t Vt, ) Versus Beam Resistance Parameter (Mp/M.) 
for Fixed Values of Pulse Duration ^ 
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Figure 38 - Plot of Duration of Initial Elastic Phase (t p/T, J Versus Pulse Duration for 

Fixed Values of Beam Resistance Parameter 
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Figure 39 - Plot of Duration of Plastic Phase Versus Beam Resistance Parameter (Mp/M^) for Fixed Values of Pulse Duration 
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Figure 40 - Plot of Duration of Plastic Phase Versus Pulse Duration (t^/T^) for Fixed Values of Beam Resistance Parameter (Mp/Mg) 
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